What constitutes a simple liquid? 
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Simple liquids are traditionally defined as many-body systems of classical particles interacting via 
radially symmetric pair potentials. We suggest that a simple liquid should be defined instead by the 
property of having strong correlation between virial and potential energy equilibrium fluctuations in 
the NVT ensemble. There is considerable overlap between the two definitions, but also some notable 
differences. For instance, in the new definition simplicity is not a property of the intermolecular 
potential because a liquid is usually only strongly correlating in part of its phase diagram. Moreover, 
according to the new definition not all simple liquids are atomic (i.e., with radially symmetric pair 
potentials) and not all atomic liquids are simple. The main part of the paper motivates the new 
definition of liquid simplicity by presenting evidence that a liquid is strongly correlating if and only 
if its intermolecular interactions may be ignored beyond the first coordination shell (FCS). This 
is demonstrated by NVT simulations of structure and dynamics of atomic and molecular model 
liquids with a shifted-forces cutoff placed at the first minimum of the radial distribution function. 
The liquids studied are: inverse power- law systems (r~" pair potentials with n = 18, 6, 4), Lennard- 
Jones (LJ) models (the standard LJ model, two generalized Kob-Andersen binary LJ mixtures, the 
Wahnstrom binary LJ mixture), the Buckingham model, the Dzugutov model, the LJ Gaussian 
model, the Gaussian core model, the Hansen-McDonald molten salt model, the Lewis- Wahnstrom 
OTP model, the asymmetric dumbbell model, and the rigid SPC/E water model. The final part of 
the paper summarizes most known properties of strongly correlating liquids, showing that these are 
simpler than liquids in general. Simple liquids as defined here may be characterized (1) chemically 
by the fact that the liquid's properties are fully determined by interactions from the molecules 
within the FCS, (2) physically by the fact that there are isomorphs in the phase diagram, i.e., curves 
along which several properties like excess entropy, structure, and dynamics, are invariant in reduced 
units, and (3) mathematically by the fact that the reduced-coordinate constant-potential energy 
hypersurfaces define a one-parameter family of compact Riemannian manifolds. No proof is given 
that the chemical characterization follows from the strong correlation property, but it is shown to be 
consistent with the existence of isomorphs in strongly correlating liquids' phase diagram. Finally, 
we note that the FCS characterization of simple liquids calls into question the basis for standard 
perturbation theory, according to which the repulsive and attractive forces play fundamentally 
different roles for the physics of liquids. 



I. INTRODUCTION 



Going back to Plato, classification or categorization is 
the epistemological process that groups objects based on 
similar properties 1\ . Having primarily biological exam- 
ples in mind, Aristotle defined categories as discrete enti- 
ties characterized by properties shared by their members 
Locke in 1690 distinguished between the nominal and 
the real essence of an object [^]. The nominal essence 
comes from experience and represents the object's ap- 
pearance, the real essence represents the object's deeper, 
constituting features. For instance, the real essence of a 
material thing is its atomic constitution, because this is 
the causal basis of all the thing's observable properties 
Q . A scientific classification is particularly useful if it re- 
flects the real essence of the objects in question by identi- 
fying their underlying common features, from which the 
more obvious and easily observable nominal properties 
follow. Having in mind Locke's concept of real essence, 
we argue below for a new definition of the class of simple 
liquids. 

Physicists love simple systems. This reflects the fun- 
damental paradigm that in order to capture a given 



phenomenon, simpler is better. Most classifications in 
physics are clear-cut, for example the classification of 
elementary particles into baryons and leptons, whereas 
classifications in other sciences usually have a wealth of 
borderline cases. Due to the diversity of molecules it is 
reasonable to expect a definition of "simple liquids" to 
be of the latter type. 

The concept of a simple liquid is old, but it remains 
central as evidenced by the 2003 book title "Basic con- 
cepts for simple and complex liquids" [5| or the review of 
non-simple liquids entitled "Theory of complicated liq- 
uids" from 2007 [6| . Generations of liquid-state theorists 
were introduced to the topic by studying Hansen and Mc- 
Donald's textbook "Theory of simple liquids" 0]. This 
book first appeared in 1976 following a period of spectac- 
ular progress in the theory of liquids, catalyzed by some 
of the very first scientific computer simulations. 

In Ref. 7 a simple liquid is defined as a classical system 
of approximately spherical, nonpolar molecules interact- 
ing via pair potentials. This and closely related defini- 
tions of liquid simplicity have been standard for many 
years ^8 12). In this definition simple liquids have much 
in common with the chemists' "nonassociated liquids" 
p^ . but there are some differences. Chemists generally 
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regard a liquid as simple even if it consists of elongated 
molecules, as long as these are without internal degrees of 
freedom and interact primarily via van der Waals forces. 
Many physicists would probably disagree. Thus it is far 
from trivial to ask: What characterizes a "simple liq- 
uid" ? More accurately: Given a classical system of rigid 
bodies with potential energy as a function of the bodies' 
centers-of-masses and their spatial orientations, is it pos- 
sible to give a quantitative criterion for how simple the 
system is? If yes, is simplicity encoded uniquely in the 
potential energy function or may the degree of simplicity 
vary throughout the phase diagram? 

Recent works identified and described the properties 
of "strongly correlating liquids" (l^ - [25j . In these liq- 
uids, by definition, the virial W and the potential energy 
U correlate strongly in their constant-volume thermal- 
equilibrium fiuctuations. Recall that the average virial 
(W) gives the contribution to pressure from intermolecu- 
lar interactions, added to the ideal-gas term NksT deriv- 
ing from momentum transport via particle motion (below 
p is the pressure, V the volume, N the number of parti- 
cles, Ub Boltzmann's constant, and T the temperature): 



pV = NkBT+{W) 



(1) 



The term "strongly correlating liquid" refers to the case 
when the WU correlation coefficient in the NVT ensem- 
ble is larger than 0.9 [l^. If angular brackets denote an 
NVT ensemble average, the correlation coefficient R is 
defined by 
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An example of a strongly correlating liquid is the stan- 
dard Lennard- Jones liquid at typical condensed- phase 
state points, i.e., not far from the solid-liquid coexis- 
tence line. Many other systems, including some molecu- 
lar models, have been shown to be strongly correlating; 
we refer the reader to papers that derive and document 
the several simple properties of strongly correlating liq- 
uids [l3 - [25j . reviewed briefiy in Ref. 25. These proper- 
ties are summarized in Sec. IIV Al after the presentation 
of the simulation results. 

The present work is motivated by developments initi- 
ated by recent findings by Berthier and Tarjus [1^ [27| . 
These authors showed that for the viscous Kob- Andersen 
binary Lennard- Jones mixture (28l . [29l | the dynamics is 
not reproduced properly by cutting the potentials at their 
minima according to the well-known Weeks-Chandler- 
Andersen (WCA) recipe [13] ■ The role of the cutoff was 
subsequently studied in two papers, showing that plac- 
ing a shifted-forces cutoff at the first minimum of the pair 
correlation functi on g ives good results for Lennard- Jones 
type systems [3l|, [sj This applies not only at moderate 
densities, but also at very high densities. A pply ing the 
same cutoff to water does not work very well [33| . Water 



is an example of a non-strongly correlating liquid with 
i? w at ambient conditions, a consequence of water's 
density maximum [l^. These findings led us to specu- 
late whether it is a general property of strongly correlat- 
ing liquids that the intermolecular interactions may be 
ignored beyond the PCS without compromising accuracy 
to any significant extent. The main part of the present 
paper shows that, indeed, using such an "FCS cutoff" 
gives accurate simulation results if and only if the liquid 
is strongly correlating. 

The paper presents results obtained from computer 
simulations of 15 different systems, some of which are 
strongly correlating. We investigate the role of the FCS 
in determining liquid structure and dynamics. Structure 
is probed by the radial distribution function (RDF), dy- 
namics by the incoherent or coherent intermediate scat- 
tering function (ISF) at the wavevector of the static 
structure factor maximum. The numerical evidence is 
clear. By varying the cutoff of the intermolecular forces, 
we find that in order to get accurate simulation results 
it is enough to take into account the interactions within 
the FCS if and only if the liquid is strongly correlating. 
In other words, for strongly correlating liquids interac- 
tions beyond the FCS are unimportant, and this applies 
only for these liquids. At present there is no rigorous ar- 
gument for this empirical "FCS property" , but we show 
that it is consistent with known properties of strongly 
correlating liquids. 

The FCS property of strongly correlating liquids docu- 
mented below emphasizes further that these are simpler 
than liquids in general. A number of other simple prop- 
erties of stron gly correlating liquids were reported pre- 
viously (l3 - l22r [2^ [25j . Altogether, these facts motivate 
our new definition of liquid simplicity. 

Section |TT] presents the results from molecular dynam- 
ics simulations and Sec. IIIII summarizes the results. Sec- 
tion IIVI gives an overview of the many simple properties 
of strongly correlating liquids, motivating our suggestion 
that a liquid is to be defined as simple whenever it is 
strongly correlating at the state point in question. Sec- 
tion |V] gives a few concluding remarks. 



II. MOLECULAR DYNAMICS SIMULATIONS 
OF ATOMIC AND MOLECULAR LIQUIDS 

In a computer simulation the intermolecular interac- 
tions, which usually extend in principle to infinity, are 
truncated at some cutoff distance beyond which they 
are ignored 3l|. To avoid a discontinuity in the force, 
which can severely affect the simulation results jiL, j5^] , 
the simulations reported below use potentials truncated 
such that the force goes continuously to zero at Tc- This 
is done by applying a so-called shifted forces (SF) cutoff 
[33 - l36j where, if the pair potential is v{r) and the pair 
force is /(r) = —v'{r), the shifted force is given by 
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/sF(r-) = 



fir) -fire) if r<rc, 







if r > Tf, 



(3) 



the case for systems with inverse power-law (IPL) pair 
potentials (w(r) oc r~"), but note that potentials with 
non-trivial angular dependence may also be Euler homo- 
geneous. 



This corresponds to using the following pair potential be- 
low vsFir) — vir)—v'{rc)ir — rc)—virc). Using an SF 
cutoff gives more accurate results and better numerical 
stability than using the standard shifted-potential (SP) 
cutoff [3l|. This is so despite the fact that an SF cutoff 
does not have the correct pair force for any r, whereas 
the pair force is correct below Tc for an SP cutoff. Ap- 
parently, avoiding discontinuity of the force at Tc is more 
important than maintaining the correct force. It was re- 
cently discussed why adding a linear term to the pair 
potential affects neither structure nor dynamics to any 
significant extent [l^. The reason is that, when one 
nearest-neighbor distance decreases, others increase in 
such a way that their sum is virtually constant. This 
argument is exact in one dimension and holds to a good 
approximation in 3D constant-volume ensemble simula- 
tions ^19i | (in constant-pressure ensemble simulations the 
volume fluctuates and the argument no longer applies). 
Coulomb interactions have also been treated by the SF 
cutoff procedure. Although the Coulomb interaction is 
long-ranged and conditionally convergent, when rc is suf- 
ficiently large an SF cutoff gives results close to those 
of the standard, much more involved Ewald summation 
method (HIHl. 

All simulations were performed in the NVT ensemble 
with periodic boundary conditions and the Nose-Hoover 
algorithm [39l - l4l| . We used the RUMD molecular dy- 
namics package developed in-house, which is optimized 
for simulations of small systems on state-of-the-art GPU 
hardware (a few thousand entities) [42j. For the molecu- 
lar models bond lengths were held fixed using the time- 
symmetrical central difference algorithm |43l-l45j. 

In the following we investigate for several systems 
whether it is possible to choose an FCS cutoff, i.e., a 
cutoff at the first minimum of the pair correlation func- 
tion, and still get the correct physics. We start by study- 
ing strongly correlating atomic liquids. Then data are 
presented for atomic liquids that are not strongly cor- 
relating. Finally data are given for two strongly corre- 
lating molecular liquids and a water model. Details of 
the models studied, number of particles, etc, are given in 
Appendix 



We simulated single-component IPL pair-potential sys- 
tems with exponents n — 18,6,4 at density p — 0.85. 
Each system was simulated at two temperatures. The 
simulated systems range from very harsh repulsive (n = 
18) to quite soft and long ranged (n = 4). Each system 
was simulated with four SF cutoffs. The role of the cut- 
off is investigated by choosing three different cutoffs: one 
placed at the first minimum of the RDF, one correspond- 
ing to the half height of the RDF from its first maximum 
to the first minimum, and one placed to the right of the 
RDF first minimum displaced the same amount as the 
difference between the first and the second cutoff. 

The RDFs g(r) are shown for n — 18, 6, 4 in Fig. [TJ 
n — 12 gives similar results (not shown). The simula- 
tions with an SF cutoff at the first minimum of the RDF 
- referred to as FCS-cutoff simulations - give a faithful 
representation of the structure. The insets show the devi- 
ations in RDF as functions of the cutoff between results 
for an SF cutoff and "true" results, quantified by inte- 
grating the numerical difference in the pair correlation 
function. Clearly, deviations increase sharply when the 
cutoff enters the FCS (blue crosses). 

We simulated also the n = 3 and n — 1 IPL fluids 
(results not shown), the latter is usually termed the one- 
component plasma. For both systems an FCS cutoff does 
not lead to the correct physics. This is most likely be- 
cause the potentials are too long ranged for such a short 
cutoff to make good sense. In fact, both are so long 
ranged that they do not have a proper thermodynamic 
limit for which the exponent must be larger than the di- 
mension [4^. An indication that an FCS cutoff works 
poorly when the IPL exponent approaches the dimen- 
sion is seen for the n = 4 simulation, for which the WU 
correlation coefficient for the FCS system starts to devi- 
ate significantly from unity. Moreover, almost invisible 
in the figure is the fact that the n — A pair correlation 
function's first maximum deviates slightly between FCS 
and true simulations. 
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A. Three inverse-power-law (IPL) fluids 

We consider first systems with 100% correlation be- 
tween virial and potential energy equilibrium fiuctuations 
in the NVT ensemble. It follows from the definition of 
the virial W = —l/3j2^i'^iU [131 that a necessary and 
sufficient condition for W to correlate perfectly with U 
is that the potential energy is a so-called Euler homo- 
geneous function of the particle coordinates r^. This is 
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FIG. 1. RDFs for single-compolient IPL fluids with exponents 
n = 18, 6, 4, each simulated at two temperatures at density 
p = 0.85. The black and orange curves show reference simu- 
lation results with large cutoffs representing the "true" IPL 
behavior, the red and green dots give results from simula- 
tions with an PCS cutoff (marked by the vertical red dashed 
lines). The insets quantify the deviations in the RDP from the 
reference RDF as functions of the cutoff; deviations increase 
dramatically when the cutoff enters the PCS (blue crosses). 
In each subfigure the virial/potential-energy correlation co- 
efficient R is given for the simulated cutoff (this quantity is 
exactly unity for IPL systems with infinite cutoff). 
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FIG. 2. Incoherent intermediatte scattering functions (ISPs) 

for the IPL fluids at the lowest-temperature state points of 

Fig. [T] The black curves give results for a large cutoff, the 

red crosses for an PCS cutoff (marked by the vertical red 

dashed lines in Fig. [T]). (a) n = 18, T = 0.30; (b) n = 6, 

T = 0.15; (c) n = 4, T = 0.10. 



B. Lennard- Jones type liquids 

Next, we consider the most studied potential in the 
history of computer simulations, the Lennard- Jones (LJ) 
pair potential. 



Figure [2] shows the incoherent intermediate scattering 
functions (ISFs) for the low-temperature state points of 
each of the three IPL systems evaluated at the wavevector 
corresponding to the first maximum of the static struc- 
ture factor. A good representation of the dynamics is 
obtained for all systems when the PCS cutofi" is used. 
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Here a and e, define, respectively, the length and energy 
scale of the interaction (dimensionless units defined by 
a = e = 1 are used below). This potential does not 
have 100% virial/potential-energy correlation, but has 
still quite strong correlations with correlation coefficients 
R > 0.9 in the condensed-fluid part of the phase dia- 
gram (also in the crystalline phase fl^). We studied 
the single-component LJ (SCLJ) liquid, two generalized 
80/20 Kob- Andersen binary LJ (KABLJ) mixtures with 
repulsive exponent 12 and attractive exponents n = 4, 10, 
and the Wahnstrom 50/50 binary LJ (WABLJ) mixture 
(Fig. [3] and Appendix [K\ gives model details). The in- 
fluence of an SF cutoff on simulation accuracy was in- 
vestigated recently for the SCLJ liquid and the standard 
KABLJ mixture {n = 6) [3l|, [s^l , but for completeness 
we include results for the SCLJ system here as well. 
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FIG. 3. The AA-particle generalized Kob- Andersen (KABLJ) 
pair potentials with fixed repulsive exponent 12 and three 
different attractive exponents 71 — 4,6, 10. The model pa- 
rameters are: eaa ~ 1-00, eas = 1.50, and ess = 0.50, 
CTAA = 1.12, aAB = 0.90, (JBB = 0.99, ruA = ms = 1. 



The role of the cutoff is again investigated by choosing 
three different cutoffs: one placed at the first minimum 
of the RDF (red color in Figs. HUH]), one corresponding 
to the half height of the RDF from maximum to the first 
minumum (blue color in Figs. HUH]), and one displaced 
to the right of the minimum the same amount as the 
difference between the first and the second cutoff (green 
color in Figs. I4ll8)) . 

In Fig. SlRDFs are shown for the SCLJ hquid at three 
different state points. The red colored circles show re- 
sults from simulations with an FCS cutoff (marked by the 
vertical red dashed line), the black curves show the cor- 
responding simulations with a large cutoff (reference sys- 
tem) . The insets quantify the deviations in the simulated 
RDF from the reference RDF as a function of the cutoff. 
The reference RDF of Figs. Sl^a) and (b) is clearly repre- 
sented well using an FCS cutoff, while choosing the cutoff 
inside the FCS results in significant deterioration. At low 
densities [Figs. IHc)], deviations occur between FCS cut- 
off simulations and the reference system. As mentioned, 
the SCLJ liquid is strongly correlating in large parts of its 
phase diagram, but as density is lowered, the correlations 
decrease gradually, and the liquid is no longer strongly 
correlating at state point (c) where R — 0.50. These 
simulations suggest that only when a liquid is strongly 
correlating is it possible to ignore interactions beyond 
the 
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FIG. 4. RDFs for the single-component Lennard- Jones 
(SCLJ) liquid at three different state points: (a) p = 0.85, 
T = 0.70 {R = 0.96); (b) p = 0.85, T = 1.00 {R = 0.97); (c) 
p = 0.55, T = 1.13 [r = 0.50). The black curves show refer- 
ence simulations with large cutoffs, the red dots/curve show 
results from simulations with an FCS cutoff (marked by the 
vertical red dashed lines) . The insets quantify the deviation in 
RDF from the reference RDF as functions of the cutoff. At all 
three state points deviations increase significantly when the 
cutoff enters the FCS (blue crosses in the insets). For state 
points (a) and (b), which are strongly correlating {R > 0.9), 
an FCS cutoff leads to accurate results. This is not the case 
for state point (c) that is not strongly correlating. 



Next, we investigated the SCLJ dynamics at the three 
state points of Fig. 2] The dynamics is studied via the 
ISF. The ISFs are shown in Fig. [5] at all state points the 
dynamics is represented well using an FCS cutoff. 
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The n = 4 case is below our definition of a strongly cor- 
relating liquid R > 0.9 (at the simulated state point this 
system has negative virial), but is still pretty well cor- 
relating. For all investigated state points an FCS cutoff 
gives accurate results. We found the same using the stan- 
dard repulsive exponent rt = 6 [33 | (results not shown). 



TABLE L Simulation results for five state points of the SCLJ 
liquid. For each state point is given density, temperature, 
correlation coefRcient, maximum deviation from the true RDF 
using an FCS cutoff, and maximum deviation from the true 
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FIG. 5. ISFs for the SCLJ liquid at the state points of Fig. 
|4l The black curves give the reference cutoff results, the red 
curves the FCS cutoff results, the blue curves results for a 
cutoff at the half-height towards the first maximum of the 
RDF, the green curves results for a cutoff to the right of the 
minimum, (a) p = 0.85, T = 0.70 {R = 0.96); (b) p = 0.85, 
T = 1.00 (7? = 0.97); (c) p = 0.55, T = 1.13 (R = 0.50). 



FIG. 6. RDFs for generalized KABLJ mixtures with repulsive 
exponent 12 and attractive exponents n = 4, 10. The black 
curves give the reference cutoff results, the red curves the 
FCS cutoff results. The insets quantify the deviation in RDF 
from the reference RDF as functions of the cutoff, (a) n = 4, 
p = 1.20, T = 0.40 {R = 0.84); (b) n = 10, p = 1.20, T = 1.00 
(ii = 0.97). 



Table [J summarizes results for five state points of the 
SCLJ liquid. The third state point is slightly below the 
triple point, i.e., in a slightly metastable state. The devi- 
ations clearly increase as the WU correlations decrease. 

We proceed to investigate mixtures of two different 
particles (A and B) interacting with LJ type potentials. 
The cutoff used for all three interactions (AA, AB, BB) 
are placed at the same distance, referring to aAA- In 
Fig. preference and FCS cutoff results are shown for the 
v4A-particle RDFs of generalized KABLJ mixtures with 
repulsive exponent 12 and attractive exponents n = 4, 10. 



The A-particle ISFs for the state points of Fig. [6] are 
shown in Fig. [71 For the KABLJ mixture, placing the 
cutoff inside the FCS (blue curves) fails to reproduce the 
dynamics properly, whereas it is well approximated using 
an FCS cutoff (red). Slight deviations are noted for the 
red curves, an issue considered in Appendix IB] that dis- 
cusses alternatives for delimiting the FCS. Similar results 
are found for the B particles (results not shown). 
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In summary, for all LJ type systems whenever there 
are strong virial potential-energy correlations (i? > 0.9), 
an FCS cutoff gives accurate results for both structure 
and dynamics. In the parts of the phase diagram where 
WU correlations are weaker, an FCS cutoff gives poor 
results. 
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FIG. 7. ISFs for generalized KABLJ mixtures with repulsive 
exponent 12 and attractive exponents n = 4, 10. The red 
and black curves give, respectively, results for FCS cutoffs 
and large reference cutoffs, (a) n = 4, p = 1.20, T = 0.40 
{R = 0.84); (b) n = 10, p = 1.20, T = 1.00 {R = 0.97). 

We also simulated the Wahnstrom 50-50 binary LJ 
mixture [47| . finding again that whenever i? > 0.9 struc- 
ture and dynamics are well reproduced using an FCS 
cutoff. We do not show these results, but show instead 
results for the coherent scattering function at one state 
point (Fig. [SJ. Again, the FCS cutoff (red crosses) gives 
the correct dynamics whereas reducing the cutoff further 
does 
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FIG. 8. Intermediate coherent scattering function for the 
Wahnstrom 50/50 binary LJ liquid at the wavevector cor- 
responding to the static structure factor maximum. The red 
and black curves give, respectively, results for a FCS cutoff 
and large reference cutoff. 



C. Buckingham liquid 



Next, we consider the single-component Buckingham 
liquid (SCB). The Buckingham potential 48, 49] is sim- 
ilar to the LJ potential, but does not have an IPL repul- 
sive term; instead the potential's short-distance behavior 
follows a steep exponential (Fig. 19]). The parameters 
of the Buckingham potential (Appendix were chosen 
such that the LJ potential is well approximated in the 
repulsive region (Fig. [3]). 
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FIG. 9. The Buckingham pair potential (red) and the LJ 
pair potential (black). The parameters of the Buckingham 
potential were chosen such that the LJ potential is well ap- 
proximated in the repulsive region. 



Figures [TOl a) and (b) show, respectively, the RDF and 
ISF for the SCB Uquid. The FCS cutoff works well. 



60 2.0 




5 





l.Oi 




0.8 










o" 


0.6 






II 






0.4 












0.2 




0.0 



1 1 






_ = 4.900 








«r^= 1.792 








xt = 1.475 (min) 


\ 






,<r^= 1.158 


\ 


R(r,= 


4.900) = 0.99 




R(r,= 


1.475) = 0.99 


Single-component Buckingham 








' (p = 1.00, T= 1.00) 










V 




^ (b) " 


1 1 






1 



10 



FIG. 10. The effect on structure and dynamics of varying 
the cutoff for the single-component Buckingham (SCB) liquid. 
The red and black curves give, respectively, results for an PCS 
cutoff and a large reference cutoff, (a) RDF at p = 1.00 and 
T = 1.00 {R = 0.99). The inset quantifies the deviation in 
RDF from the reference RDF as a function of the cutoff, (b) 
ISF at the same state point. 
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D. Dzugutov liquid 
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Figure [TT] shows the Dzugutov (DZ) pair potential [50| . 
which was originally suggested as a model potential im- 
peding crystallization by energetically punishing particle 
separations corresponding to the next-nearest neighbor 
distance of a face-centered cubic lattice. At short dis- 
tances the DZ pair potential approximates the LJ poten- 
tial. 
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FIG. 12. The effect on structure and dynamics of varying the 
cutoff for the Dzugutov (DZ) liquid at p = 0.80 and T = 0.75 
{R — 0.71). The red and black curves give, respectively, re- 
sults for an FCS cutoff and a large reference cutoff, (a) RDF 
where the inset quantifies the deviation from the reference 
RDF as a function of the cutoff, (b) The intermediate coher- 
ent scattering function at the same state point, including here 
results for a cutoff within the FCS (blue crosses). 



FIG. 11. The Dzugutov (DZ) pair potential (blue). Also 
shown is the LJ pair potential (black curve). The DZ potential 
approximates the LJ potential around the first minimum, but 
has a local maximum at larger distances. 



E. Lennard-Jones Gaussian liquid 



Figures fT2la) and (b) show, respectively, the RDF and 
intermediate coherent scattering function of the DZ sys- 
tem. For this system the use of an FCS cutoff leads to 
poor results. This is not surprising given the fact that 
using an FCS cutoff removes the maximum of the DZ 
potential. What is important here, however, is that the 
poor FCS cutoff results correlate with the fairly weak 
virial/potential-energy correlations {R = 0.71). 



The Lennard-Jones Gaussian liquid [51| (LJG) is a 
non-strongly correlating liquid with the two-minimum 
pair potential shown in Fig. [131 The parameters of LJG 
model (Appendix |A| are such that the second LJG po- 
tential minimum does not coincide with that of the SCL J 
system [52|. 
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FIG. 13. The Lennard-Jones Gaussian (LJG) pair potential 
[sH constructed by adding a Gaussian to the pair LJ poten- 
tial. Two distinct minima are present. 



Results from simulating structure and dynamics of the 
LJG hquid are shown in Figs. [HUa) and (b). The FCS 
cutoff does not give the correct RDF, whereas deviations 
in the dynamics are fairly small. Note that the FCS cutoff 
removes the second minimum. 
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F. Gaussian core model 

The Gaussian core model (GCM) [H [111, which is 
not strongly correlating, is defined by a Gaussian pair 
potential and thus has a finite potential energy at zero 
separation. The high-density regime of the GCM model 
(p > 1.5) has recently received attention as a single- 
component model glass former [ssl ] because it is not prone 
to crystallization and shows the characteristic features of 
glass-forming liquids (large viscosity, two-step relaxation, 
etc). 

Figure US] shows the RDF and ISF for the GCM liquid. 
The GCM crystallizes when an FCS cutoff is used. For 
this reason, obviously, an FCS cutoff is not able to repro- 
duc 
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FIG. 15. The effect on structure and dynamics of varying the 
cutoff for the GCM liquid. The red and black curves give, 
respectively, results for an FCS cutoff and a large reference 
cutoff, (a) RDF at p = 0.85 and T = 0.01 {R = -0.77). 
The inset quantifies the deviation in RDF from the reference 
RDF as a function of the cutoff. The red curve represents a 
crystallized state, (b) ISF at the same state points. 



G. The Hansen-McDonald molten salt model 



FIG. 14. The effect on structure and dynamics of the cutoff for 
the LJG liquid. The red and black curves give, respectively, 
results for an FCS cutoff and a large reference cutoff, (a) RDF 
at p = 0.85 and T = 2.00 {R = 0.44). The inset quantifies 
the deviation in RDF from the reference RDF (black curve) 
as a function of the cutoff, (b) ISF at the same state point. 



The final atomic system studied is the "singly-charged 
molten salt model" of Hansen and McDonald [5^. In 
Fig. [16] we see that the structure is not represented well 
by the use of an FCS cutoff. Interestingly, the dynamics 
is well reproduced using this cutoff - even better, in fact, 
than for a larger cutoff (Fig. [TCTb)). 
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FIG. 16. The effect on structure and dynamics of varying 
the cutoff for Ifansen-McDonald singly-charged molten salt 
model. Because of the competing interactions (Coulomb and 
n — 9 repulsive IPL) this model is not strongly correlating. 
The red and black curves give, respectively, results for an PCS 
cutoff and a large reference cutoff, (a) AA-particle RDF at 
p = 0.37 and T = 0.018 {R = 0.15). The inset quantifies the 
deviation in RDF from the reference RDF as a function of the 
cutoff, (b) y4-particle ISF at the same state point. 



FIG. 17. The effect on structure and dynamics of varying the 
cutoff for the Wahnstrom OTP model. The red and black 
curves give, respectively, results for an FCS cutoff and a large 
reference cutoff, (a) RDF of the LJ particles at p = 0.33 and 
T — 0.70 {R = 0.91). The inset quantifies the deviation in 
RDF from the reference RDF as a function of the cutoff. The 
spikes derive from the bonds, (b) ISF at the same state point. 



H. Two strongly correlating molecular model 
liquids 



We finish the presentation of the numerical results by 
giving data for three molecular model liquids. In this sub- 
section data are given for two strongly correlating molec- 
ular liquid models, the Lewis-Wahnstrom OTP [s^, 
and the asymmetric dumbbell [24] models, which rep- 
resent a molecule by three and two rigidly bonded LJ 
spheres, respectively. The next subsection gives data for 
a rigid water model. 

Figures [TTT a) and (b) show the LJ particle RDF and 
ISF of the OTP model. Both quantities are well approx- 
imated using an FCS cutoff, although slight deviations 
are noted for the ISF (red curve, see Appendix B for 
considerations concerning this). The OTP model is a 
border-line strongly correlating liquid (R = 0.91). 



Figures ITSTa) and (b) show corresponding figures for 
the large (A) particle of the asymmetric dumbbell model 
at a viscous state point. The use of an FCS cutoff gives 
accurate results for both structure and dynamics. The 
FCS cutoff was placed at the second minimum of the AA 
RDF, because the AA RDF has here a lower value than 
at the first minimum. If the cutoff is placed at the first 
minimum, clear deviations are noted (data not shown). 
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FIG. 18. The effect on structure and dynamics of varying 
the cutoff for the asymmetric dumbbell model. The red and 
black curves give, respectively, results for an PCS cutoff and 
a large reference cutoff. Note that in this case the FCS cutoff 
is defined by using the second minimum (the first minimum 
is not the absolute minimum), (a) RDF at p = 0.93 and 
T = 0.46 {R = 0.96). The inset quantifies the deviation in 
RDF from the reference RDF (black curve) as a function of 
the cutoff, (b) A-particle ISF at the same state point. 



I. Rigid SPC/E water 



We consider finally the rigid SPC/E water model [s^ 
(Fig. \T9\\ . This model is not strongly correlating at am- 
bient conditions, a fact that reflects water's well-known 
density maximum 19|. The structure of the SPC/E wa- 
ter model is not well represented using an FCS cutoff. 
The FCS cutoff dynamics on the other hand shows only 
slight deviations from that of the reference curve (black). 
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FIG. 19. The effect on structure and dynamics of varying the 
cutoff for the rigid SPC/E water model ^59]. The red and 
black curves give, respectively, results for an FCS and a large 
reference cutoff, (a) Oxygen-oxygen RDF at p = 1.00 and 
T = 4.00 {R = 0.08). The inset quantifies the deviation in 
RDF from the reference RDF as a function of the cutoff, (b) 
Oxygen ISF at the same state point. 



III. SUMMARIZING THE SIMULATION 
RESULTS 

The previous section showed that structure and dy- 
namics are well approximated in simulations using an 
FCS cutoff for the following atomic and molecular sys- 
tems: 

• Inverse power-law (IPL) systems (n — 18,6,4), 

• single-component Lennard- Jones liquid at density 
p = 0.85, 

• generalized Kob-Andersen binary Lennard- Jones 
mixtures, 

• Wahnstrom binary Lennard- Jones mixture, 

• single-component Buckingham liquid, 

• Lewis- Wahnstrom OTP model, 

• asymmetric dumbbell model. 

These systems are all strongly correlating p^| - |22| . Thus 
for strongly correlating liquids it is enough to know the 
intermolecular interactions within the FCS in order to 
accurately simulate structure and dynamics. 

The simulations showed further that for all of the fol- 
lowing atomic and molecular systems structure and/or 
dynamics are not properly reproduced when an FCS cut- 
off is used: 

• Dzugutov (DZ) liquid, 

• Lennard- Jones Gaussian (LJG) liquid, 

• Gaussian core model (GCM), 

• Hansen-McDonald molten salt model. 
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• rigid SPC/E water model. 

None of these liquids are strongly correlating. For all 
these systems larger cutoffs are needed in order to faith- 
fully reproduce the system's physics. 

In conclusion, a shifted-forces FCS cutoff leads to accu- 
rate results if and only if the liquid is strongly correlating 
at the state point in question. We know of no exceptions 
to this rule. This suggests that strongly correlating liq- 
uids are characterized by the property that intermolecu- 
lar interactions beyond the FCS can be ignored. 



IV. THE REAL ESSENCE OF LIQUID 
SIMPLICITY 



As mentioned in the introduction a definition of sim- 
le liquids is most useful if it identifies their real essence 
, the underlying fundamental characteristic from which 
these liquids' simple features, their nominal essences, fol- 
low. We suggest below that the class of simple liquids is 
to be identified with the class of strongly correlating liq- 
uids (Sec. IIV Cp . This is motivated by first summarizing 
the many simple properties of strongly correlating liq- 
uids (Sec. IIV Ap . then showing that this class of liquids 
can be characterized from three different perspectives: 
mathematically, physically, and chemically (Sec. IIV Bp . 
This gives three very different but equivalent character- 
izations, indicating that the class of strongly correlat- 
ing liquids is fundamental and further motivating the 
suggestion that the real essence of liquid simplicity is 
the existence of strong correlations of virial/potential- 
energy equilibrium NVT fluctuations. By connecting to 
the chemists' concept of non-associated liquids we then 
discuss which real- world liquids are simple (Sec. lIVPp . 
discuss briefly liquids near interfaces (Sec. lIVEp . and 
give examples of complex liquid properties (Sec. lIVFp . 
Finally, Sec. IIVGI points out that our results call into 
question traditional perturbation theory, which is based 
on assuming quite different roles of the attractive and the 
repulsive forces. 



A. Strongly correlating liquids and their simple 
properties 

Most properties of strongly correlating liquids follow 
from the existence of "isomorphs" in their phase diagram 
(see below). Some of their properties were identified be- 
fore isomorphs were defined in 2009 [2l|, however, for 
instance that 

• all eight fundamental thermoviscoelastic response 
functions are given in terms of just one, i.e., the 
dynamic Prigogine-Defay ratio is close to unity [l7| , 

• aging may be described by adding merely one extra 
parameter (T9l.[20j. 



power-law density scaling [60] is obeyed to a good 
approximation, i.e., for varying density and tem- 
perature the relaxation time is a function of /T 



2A\. 



An isomorph is an equivalence class of state points in 
a system's phase diagram. Two state points (pi, Ti) 
and (/92, T2) are defined to be isomorphic if the follow- 



ing holds [21| : Whenever one microconfiguration of state 
point (1) and one of state point (2) have the same re- 
duced coordinates (i.e., p^'^r'^'' = ^2^'^''^'' fo'" parti- 
cles i), these two microconfigurations have proportional 
configurational Boltzmann factors. 



(5) 



For most systems the isomorph concept is approximate 
just as WU correlations are rarely perfect. Thus we do 
not require Eq. ([5]) to be rigorously obeyed for all mi- 
croconfigurations, but only to a good approximation and 
only for all physically relevant microconfigurations. By 
this is meant microconfigurations that are not a priori 
unimportant for the physics. Being an equivalence class, 
an isomorph defines a continuous curve of state points in 
the liquid's phase diagram. Only liquids for which the 
potential energy is an Euler homogeneous function, e.g., 
systems with IPL pair potentials, have exact isomorphs. 
An IPL fluid with v{r) cx r~" has exact isomorphs char- 
acterized hj /T = Const, where 'y = nl2>. 

Appendix A of Ref. [2l| showed that a liquid is strongly 
correlating if and only if it has isomorphs to a g ood ap- 
proximation. This was confirmed in Refs. [2l| and |22| . 
which showed that Lennard- Jones type atomic liquids 
have good isomorphs. Reference [6l| showed that the 
strongly correlating Lewis- Wahnstrom OTP and asym- 
metric dumbbell molecular models also have good iso- 
morphs. 

Equation ([5]) has many consequences. These were de- 
rived and discussed in detail in the original isomorph 
paper from 2009 (Ref. i21,), to which the reader is re- 
ferred. Basically, structure and dynamics at two isomor- 
phic state points are identical in reduced units. Quan- 
tities that are invariant along an isomorph include (but 
are not limited to): 

1. The excess entropy, i.e., the entropy in excess of 
the ideal gas entropy at the same density and tem- 
perature - this is the configurational contribution 
to the entropy (a term that is negative because a 
liquid is always more ordered than an ideal gas at 
same density and temperature). 

2. All A^-body entropy terms. Recall that the excess 
entropy can be expanded in a series of two-body, 
three-body, etc, terms. Each term is invariant along 
an isomorph [21]. 

3. The isochoric heat capacity. 
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4. The structure in reduced units (defined by = 
p^l^Vi for all particles i). Not only the radial dis- 
tribution function, but all higher-order distribution 
functions in reduced units are isomorph invariant. 

5. The Newtonian NVE znA Nose-Hoover iVFT equa- 
tions of motion in reduced units; likewise Brownian 
dynamics. 

6. All autocorrelation functions in reduced units. 

7. All average relaxation times in reduced units. 

8. Reduced transport coefficients like the diffusion co- 
efficient, viscosity, etc. 

Isomorphs have the further interesting property that 
there is no relaxation for an instantaneous change of tem- 
perature and density when jumping from an equilibrated 
state point to a different state point isomorphic with the 
initial state. The absence of relaxation derives from the 
fact that the Boltzmann probabilities of scaled microcon- 
figurations are identical. Such "isomorph jumps" work 
quite well for the KABLJ liquid [Hj, for the asymmetric 
dumbbell, and for the Lewis- Wahnstrom OTP molecu- 
lar models [6l|. Moreover, the effective temperature of 
a glass prepared by a temperature-density jump from an 
equilibrium state of a strongly correlating liquid depends 
only on the final density [16]; this provides yet another 
example of a simple feature of these liquids. 

Some further predictions for the class of strongly cor- 
relating liquids deriving from the existence of isomorphs 
are: 

• The solid-liquid coexistence curve is an isomorph 
[21I [2^ . This implies invariance along the coex- 
istence curve of the reduced structure factor, the 
reduced viscosity, the reduced diffusion constant, 
etc, as well as pressure invariance of the melting en- 
tropy, and a reduced-unit Lindemann melting cri- 
terion [U. 

• Collapse of the two-order-parameter maps of 
Debenedetti et al. j62| - [66{ to one-dimensional 
curves [2l| . 

• Isochronal superposition [63|, i.e., the fact that 
when pressure and temperature are varied, the av- 
erage relaxation time determines the entire relax- 
ation spectrum [2l| . 

The above listed properties of strongly correlating liq- 
uids all reflect simple features of strongly correlating liq- 
uids. A final, recently established simple property is a 
thermodynamic separation identity: for all strongly cor- 
relating liquids, if s is the excess entropy per particle, the 
temperature as a function of s and density p factorizes 
as follows 



T = f{s)h{p). 



(6) 



Equation ([6]) has a number of consequences [11], includ- 
ing the Gruneisen equation of state and that the iso- 
morphs of LJ liquids - in particular, the LJ solid-liquid 
coexistence curve - are given by {Ap"^ — Bp^)/T = Const. 

mm- 



B. Mathematical, physical, and chemical 
characterizations of strongly correlating liquids 

At a given state point, if the average potential energy is 
denoted by (U) , the constant-potential-energy hypersur- 
face is defined by fi = {(ri, rjv) G R'^^\U{ri, rjv) — 
([/)}. This is a compact, Riemannian (37V — 1)- 
dimensional differentiable manifold. Each state point has 
its own such hypersurface. In this way, a family of man- 
ifolds is defined throughout the phase diagram. In Ap- 
pendix A of Ref. [21I it was shown that the reduced-unit 
constant-potential-energy manifold is invariant along a 
strongly correlating liquid's isomorphs, and that invari- 
ance curves exist for these manifolds only for strongly 
correlating liquids. Thus, for such liquids these mani- 
folds constitute a one-parameter family of manifolds, not 
two-parameter families as expected from the fact that the 
phase diagram is two-dimensional. 

The above provides a mathematical characterization of 
the class of strongly correlating liquids. The physical 
characterization was discussed already: the existence of 
isomorphs in the phase diagram. A liquid is strongly 
correlating if and only if it has isomorphs to a good ap- 
proximation (Appendix A of Ref. 21). 

The chemical characterization of strongly correlating 
liquids is the property documented in the present paper: 
A liquid is strongly correlating at a given state point if 
and only if the liquid's structure and dynamics are accu- 
rately calculated by shifted-forces cutoff simulations that 
ignore interactions beyond the first coordination shell. 
This is an empirical finding for which we have at present 
no compelling arguments, but can at present merely point 
out two things. First, the property of insignificance of in- 
teractions beyond the FCS is an isomorph invariant: //a 
liquid has good isomorphs and if an FCS cutoff works well 
at one state point, FCS cutoffs must work well for all its 
isomorphic state points. Thus the chemical characteriza- 
tion of strongly correlating liquids is consistent with the 
fact that these liquids have isomorphs. Secondly, almost 
all of the fiuctuations in virial and potential energy of 
the LJ liquid come from interparticle separations within 
the FCS [l|. 

The new "FCS characterization" of strongly correlat- 
ing liquids shows that these liquids are characterized by 
having a well-defined FCS. Most likely it is the existence 
of a well-defined FCS which implies the almost cancel- 
lation of the linear term of the shifted-force potential. 
The fact that interactions beyond the FCS may be ig- 
nored shows that interactions are effectively short ranged, 
which means that the structure is dominated by what 
may be termed packing effects. 
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C. Defining tlie class of "simple liquids" 

Section IIV Al listed strongly correlating liquids' many 
simple properties. Section IIVBI showed that this liq- 
uid class may be characterized from three quite different 
points of view. Clearly the class of strongly correlating 
liquids is fundamental. Since the properties of strongly 
correlating liquids are generally simpler than those of liq- 
uids in general, we now propose the definition 

• Simple liquids = Strongly correlating liquids 

This is the basic message of the present paper, which 
implies a quantification of the degree of simplicity via 
the number R of Eq. ([2]) , the NVT ensemble equilibrium 
virial/potential-energy correlation coefficient. 

Compared to the standard definition of simple liquids 
(as those with radially symmetric pair interactions) there 
are some notable differences: 

1. Simplicity is quantified by a continuously variable, 
it is not an on/off property. 

2. The degree of simplicity generally varies through- 
out the phase diagram. Most strongly correlating 
liquids loose this property as the critical point is 
approached; on the other hand many or most non- 
strongly correlating liquids are expected to become 
strongly correlating at very high pressures [7l|. 

3. Not all "atomic" liquids (i.e., with radially sym- 
metric pair interactions) have simple regions in the 
low-pressure part of the phase diagram (compare 
the Dzugutov, Lennard-Jones Gaussian, Gaussian 
core, and molten salt models). 

4. Not all simple liquids are atomic (compare the 
Wahnstrom OTP and the asymmetric dumbbell 
models). 

According to the new definition of liquid simplicity the 
case where the potential energy is an Euler homogeneous 
function of the particle positions [R = \) sets the gold 
standard for simplicity. This is consistent with the many 
simple properties of IPL liquids. Due to the absence of 
attractions, IPL fluids have no liquid-gas phase transi- 
tion. In this sense it may seem strange to claim that IPL 
fluids are the simplest liquids. However, more realistic 
strongly correlating liquids like the LJ liquid cease to be 
so when the liquid-vapor coexistence line is approached 
near the critical point, showing that this phase transi- 
tion cannot be understood in the framework of simple 
liquids. This contrasts with the liquid-solid phase transi- 
tion, where for instance the fact that the coexistence line 
for simple liquids is an isomorph ~ confirmed by simula- 
tions of the LJ li quid 1221 - explains several previously 
noted regularities |2l|, [63 . 

Is the hard-sphere fluid simple? One may define a 
configurational virial function for this system, but it is 
not obvious how to define a potential energy function 



that is different from zero. Thus there is no meaningful 
correlation coefficient R for hard-sphere fluids. On the 
other hand, the hard-sphere liquid may be regarded as 
the n — ^ 00 limit of an IPL liquid, and it is well known 
that for instance the hard-sphere radial distribution func- 
tion is close to that of, e.g., an r^^^ IPL liquid at a 
suitably chosen temperature. This would indicate that 
hard-sphere liquids are simple, consistent with the pre- 
vailing point of view. Another interesting case is that of 
the Weeks-Chandler- Andersen (WCA) version of the LJ 
liquid, which cuts off all attractions by putting the force 
equal to zero beyond the potential energy minimum. This 
liquid is strongly correlating [T^ , but in simulations that 
the WCALJ liquid has somewhat poorer isomorphs than 
the LJ liquid. 

It is possible that the hard-sphere liquid and the 
WCALJ liquid should be both excluded from the class 
of simple liquids on the grounds that their potentials are 
not analytic. For systems interacting with pair poten- 
tials, for instance, it could make good sense to add the 
extra requirement that the pair potential is an analytical 
function of the inverse pair distance, i.e., that an expan- 
sion exists of the form v{r) = 'Ylin'^riT~"' ■ Such an extra 
analyticity requirement would not exclude any strongly 
correlating liquids occurring in nature where all poten- 
tials are expected to be analytic. 

D. Which liquids in the real world are simple? 

Real-world liquids may be classified according to the 
nature of the chemical bonds between the molecules. 
There are five types of bonds [zH, listed below with a 
few typical examples (polymeric systems may be added 
as a separate class): 

• Van der Waals bonds (argon, toluene, butane, ...); 

• Metallic bonds (gold, aluminum, alloys, ...); 

• Hydrogen bonds (water, glycerol, ethanol, ...); 

• Ionic bonds (molten sodium chloride and potassium 
nitrate, room-temperature ionic liquids, ...); 

• Covalent bonds (silica and borate melts, ...). 

Most liquids involve elements of more than one type 
of chemical bonds. For instance, van der Waals forces 
are present in all liquids; the first class consists merely 
of those liquids that only have van der Waals forces. 
Another borderline example is a dipolar organic liquid 
like di-butyl-phthalate, where van der Waals as well as 
Coulomb forces are present; a liquid like glycerol has also 
strong dipolar interactions, i.e., an element of the ionic 
bonds that defines class 4), etc. 

Based on computer simulations and known properties 
of liquids we believe that most or all van der Waals and 
metallic liquids are strongly correlating [l3, [U m, i.e., 
simple. Liquids that are not simple are the hydrogen-. 
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ionically, and covalently bonding liquids. In these cases 
the virial/potential-energy correlations are weakened by 
one or the other form of competing interactions. 

Metals play a special role as simple liquids, because 
their interatomic forces derive from collective interactions 
between ion cores and free electrons [!]■ The resulting 
interaction is a non-directional interaction between sym- 
metric ion cores, i.e., these systems are "simple" in the 
traditional sense. Preliminary computer simulations of 
ours show that metals are strongly correlating [l^, so 
metals are simple also in the sense of this paper, as well. 
However, not all isomorph invariants apply for metals. 
For instance, the electron gas can influence the collective 
dynamics without any visible structural and relaxational 
counterpart [73 . [tsI ] , so isomorph invariance most likely 
breaks down for these (fast) collective degrees of freedom. 

It should be emphasized that the above considerations 
refer to ambient or moderate pressure conditions. It was 
recently suggested that if crystallization is avoided, all 
liquids become strongly correlating at sufficiently high 
pressure [7l|. Thus, e.g., the molten silicates of the 
earth's upper mantle are predicted to be simpler than 
molten silicates at ambient pressure. 



E. Liquids near interfaces 



F. A note on complex liquid behavior 

Liquids that are not simple in the above defined sense 
of the term often have complex properties [H, |82| - [85| . 
Water with its correlation coefficient close to zero at am- 
bient conditions is a prime example of a complex liq- 
uid. It is well known for water that a certain region of 
state points in the density /temperature phase diagram 
exhibits anomalous thermodynamic behavior in the sense 
that isobaric heating implies densification. Numerical ev- 
idence indicates that these state points lie within a larger 
region with diffusion anomaly, i.e., an increased diffu- 
sivity upon isothermal compression |64||, a region that 
in turn lies within a larger region of structural anomaly 
characterized by decreasing order upon isothermal com- 
pression (63 |. 

Different order parameters exist for characterizing the 
structural order of liquids, some of which relate purely to 
an integral over the RDF [H-dJ]. In this way it is pos- 
sible to calculate the contribution to structural anoma- 
lies from the different coordination shells (82l - [8^ . It has 
been shown [83, 84] that the structural anomaly of wa- 
ter and water like liquids is not a: "... first-shell effect. 
Rather, they reflect how structuring in second and more 
distant coordination shells responds to changes in ther- 
modynamic or system parameters." The full anomalous 
behavior of water derives from interactions beyond the 
FCS [13, [13] • This is consistent with the results pre- 
sented in this paper, since the structure and dynamics of 
strongly correlating liquids are given exclusively by the 
interactions within the FCS. 



Since the property of being simple cannot be read of 
from knowledge of the potential alone, it is of interest to 
consider liquids under more general circumstances, for 
instance under conflnement or generally near interfaces. 
Liquids near interfaces show rich and complicated be- 
havior. For instance, a liquid confined to the nanoscale 
may change its dynamic properties several orders of mag- 
nitude compared to the bulk system. Predicting these 
changes is an important challenge relevant for biological 
systems, engineered devices, etc. Recently, it was shown 
that some liquids retain bulk liquid behavior in confine- 
ment [76l - [79j . More specifically, it was shown that Rosen- 
feld's excess entropy scaling in the bulk persists in con- 
finement and is to a good approximation independent of 
the wall-fluid interaction strength and the degree of con- 
finement. This was shown for LJ and hard-sphere liq- 
uids, suggesting the possibility of extending the concept 
of a simple liquid to apply beyond bulk systems. Given 
that interactions in strongly correlating liquids are lim- 
ited in range by the radius of the FCS, one may speculate 
that these liquids are simple also by having the property 
that an external field at one point affects the liquid over 
shorter distances than for a general liquid. Clearly, more 
work is needed to clarify the relevance and consequences 
of the present definition of liq uid simplicity near inter- 
faces and in external fields [80, SjJ • 



G. To which extent do the assumptions of standard 
pertubation theory hold? 

The finding that the FCS plays a crucial role for a large 
class of systems may be taken as a modern demonstration 
of the classic van der Waals picture of liquids in the sense 
that such liquids can be understood in terms of packing 
effects ^8^. On the other hand, our results call into ques- 
tion the basis of traditional perturbation theory, which is 
also usually traced back to van der Waals [87[ . Perturba- 
tion theory is based on the assumption of entirely differ- 
ent roles being played by the repulsive and the attractive 
forces [3, [13, [IS HBUll : The repulsive forces define the 
structure and reduce the entropy compared to that of an 
ideal gas at same density and temperature, the attractive 
forces reduce the pressure and energy compared to that 
of an ideal gas. From the findings of the present and a 
previous paper [32] it is clear, however, that this picture 
applies only when the FCS coincides roughly with the 
pair potential minimum, i.e., at fairly low pressure. At 
very high pressures the whole FCS is within the range 
of the repulsive forces, here the attractive forces play lit- 
tle role for simple liquids. In general, what is important 
for a strongly correlating liquid is to take into account 
properly all forces from particles within the FCS - and 
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only these. Thus the well-known WCA reference system, 
which ignores the attractions, is a good reference only at 
such high pressures that all forces from particles within 
the FCS are repulsive [Mllllll- 

The dominance of the FCS for simple liquids reflects 
the fundamental physics that the characteristic length 
defining the potential minimum (e.g., a of the LJ poten- 
tial) is much less important than generally believed: a de- 
termines the density of the low-pressure condensed phase, 
but that is all. The only physically important length for 
simple liquids is that given by the macroscopic density 
itself: p^^/"^. At low pressures this length is roughly that 
of the potential energy minimum, explaining why the lat- 
ter has been generally assumed to be important. 

The above considerations apply only for simple liq- 
uids; in general both lengths play important roles for the 
physics. The irrelevance of any length defined by the mi- 
croscopic potential emphasizes once again that the class 
of strongly correlating liquids is at the one end of the 
"complexity scale" where, at the other end, one finds 
systems like macromolecules, electrolytes, interfaces, mi- 
celles, or enzymes, for which multiple length and time 
scales are important [ooj . 

V. CONCLUDING REMARKS 

If you ask a chemist what is a simple liquid, he 
or she may answer that nonassociated liquids are sim- 
ple, whereas associated liquids are much more complex. 
These two concepts are defined as follows in Chandler's 
textbook (isj . The intermolecular structure of a nonas- 
sociated liquid "can be understood in terms of packing. 
There are no highly specific interactions in these sys- 
tems." In contrast, water is an example of an associ- 
ated liquid, and its "linear hydrogen bonding tends to 
produce a local tetrahedral ordering that is distinct from 
what would be predicted by only considering the size and 
shape of the molecule" [T^. 

Packing usually refers to purely entropic, hard-sphere 
like behavior. Given that no realistic potentials are in- 
finitely repulsive it makes good sense to interpret packing 
more generally as all short-ranged effects of the inter- 
molecular interactions. If one accepts this more general 
interpretation of packing, the crucial role of the FCS for 
strongly correlating liquids is consistent with the under- 
standing since long that the properties of nonassociated 
liquids can be interpreted in terms of packing: 

• Once the forces from particles within the FCS are 
known, basically everything is known. 

In other words, for a simple liquid there are no impor- 
tant long-range interactions, and "considering the size 
and shape of the molecule" is enough to account for the 
liquid's physical properties. This applies even for the r~'^ 
IPL fiuid, which has a fairly long-ranged interaction. 

The present definition of the class of simple liquids 
is thus consistent with the chemists' general picture of 



simple liquids. The new definition goes further, how- 
ever, by quantifying simplicity via the virial/potential- 
energy correlation coefficient R of Eq. In particu- 

lar, the degree of simplicity is not an on/off property of 
the potential, but varies continuously with state point. 
Thus even a complex liquid like water is expected to ap- 
proach simple behavior under sufficiently high pressure 
[Til and, conversely, the prototype strongly correlating 
LJ liquid becomes gradually more complex as density 
is lowered and the critical region and/or the gas phase 
is approached. Is this a problem, given that everyone 
agrees that the gas phase is simple? We do not think 
so. In fact, the gas phase is simple for a different rea- 
son, namely that molecules move freely most of the time, 
only interrupted by occasional fast and violent collisions 
with other molecules. It would be strange if a system 
exhibiting one form of simplicity could be transformed 
continuously, while maintaining its simplicity, into a sys- 
tem of an entirely different form of simplicity; one would 
expect the intermediate phase to be complicated. 

Liquid simplicity is characterized by the correlation 
coefficient R of Eq. 1^ being close to unity, i.e., that 
1 — i? is a small number. This situation is typical in 
physics, where simplifying features always appear when 
some dimensionless number is small. The obvious ques- 
tion arises whether a perturbation theory may be con- 
structed around simple liquids, embracing the more com- 
plex ones. Only future can tell whether this is possible, 
but it does present a challenge because the properties of 
IPL fiuids {R — 1) cannot be worked out analytically. 

A potentially annoying feature about defining liquid 
simplicity from the existence of strong correlations of the 
virial/potential-energy fluctuations is that one cannot 
read off directly from potential and state point whether 
or not a given liquid is simple. The same holds for the ex- 
istence of isomorphs or the property that an FCS cutoff 
reproduces the correct physics. We believe one should 
accept this as an acceptable cost for precisely delimit- 
ing the simple liquids from others. With the power of 
today's computers this is much less of a problem than 
previously. For most systems a brief simulation will de- 
termine whether or not the liquid is strongly correlating 
at the state point in question. 

Except for the IPL fluids, no system is simple in the 
entire fluid phase. This paper focused on the condensed 
liquid phase, not too far from the solid-liquid coexistence 
line, but far from the critical point and the gas phase - 
it is here that some liquids are simple. The focus on 
liquids is not meant to imply a limitation to the liquid 
phase, however. In fact, simulations show that when a 
strongly correlating liquid crystallizes, the crystal is at 
least as strongly correlating A theory has been de- 
veloped for (classical) strongly correlating crystals, show- 
ing that the property of strong virial/potential-energy 
equilibrium fluctuations in the NVT ensemble is an an- 
harmonic effect that survives as T — )■ [l^. Of course, 
low-temperature crystals are not classical systems, and 
both for liquids and crystals an interesting topic for fu- 
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ture works is the implication of the proposed simphcity 
definition for the quantum description of condensed mat- 
ter. 

Section [IV Al summarized the several nominal essences 
of simple liquids. What is the real essence of liquid 
simplicity? Given that three fundamental characteriza- 
tions of strongly correlating liquids are equivalent - the 
mathematical, the physical, and the new chemical (FCS) 
characterization - this question cannot be answered un- 
equivocally. At the end of the day it is a matter of 
taste whether one defines liquid simplicity from the exis- 
tence of strong virial/potential-energy correlations, from 
the existence of isomorphs, from the existence of invari- 
ance curves for constant-potential-energy hypersurfaces, 
or from the property that interactions beyond the FCS 
play little role. All four properties are equivalent. 
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Appendix A: Model details 

The model systems investigated are listed below. 
Quantities are given in rationalized units defined by set- 
ting e = cr = 1. Masses that are not specified are unity. 
Single- component inverse power-law (IPL) fluids: N = 
1024 particles interacting via v(r) = e(cr/r)". Three dif- 
ferent fluids were studied (n = 18, 6, 4). 
Single- component Lennard-Jones liquid: N = 1024 par- 
ticles interacting via Eq. 

Generalized Koh- Andersen binary mixture [28l[29|: A bi- 
nary mixture of 820 A particles and 204 B particles in- 
teracting via 

v{r) = e„^/(12 - n) [nia^^p/rf^ - 12K^/r)"] . Bi- 
nary mixtures with n = 4, 10 were studied. The param- 
eters used are: e^A = 1.00, eab = 1-50, cbb = 0.50, 
aAA = 21/6, aAB = 0.8 • 2"/^ ubb = 0.88 • 2^1^. 
Buckingham liquid: N = 1000 

particles interacting via v{r) — 
e [6/(a - 6) exp[a(l - r/r,n)] -a/{a- 6)(r„/r)'^]. 
The parameters used are 6=1,0; = 14.5, = 2^/^. 
Dzugutov liquid [sot : N — 1024 particles interacting via 
v{r) ^ vi + V2 where vi — (^(r^" — _B)exp(c/(r — a))) 
and V2 = B exp{d/{r — b)) and r>a^vi — 0, r> 
6 ^ W2 = (a < 6). The parameters used are a = 1.87, 
b = 1.94, c = 1.1, d = 0.27, A = 5.82, B = 1.28, n = 16. 
Lennard-Jones Gaussian liquid [sij : N ~ 
1024 particles interacting via v(r) = 
e (((T/r)i2 - 2{a/rf + eo exp[-(r - ro)V2a2]) The 
parameters used are ctq = 0.02, eo = 1.50, rg = 1.47. 



Gaussian core model [5l'|: N = 1024 particles interacting 
via v{r) = eexp [— (t'/ct)^] . 

The Hansen-McDonald molten salt model [1^: A'' = 
2744 particles forming an equimolar binary mixture of 
singly charged cations and anions. The potential be- 
tween two particles of charge qa and qp is given by 
v{r) = (l/9)r~^ + qaqp/r, where q+ ^ q_ ~ -1. 

Lewis- Wahnstrom OTP [H^jlHl]: This model consists of 
three identical LJ particles rigidly bonded in an isosceles 
triangle with unity sides and top-angle of 75° (number of 
molecules studied: N — 320). 

The asymmetric dumbbell model .24]: This model con- 
sists of a large (A) and a small (B) LJ particle, rigidly 
bonded with bond distance r^s = 0.29/0.4963 (number 
of molecules studied: 500). The asymmetric dumbbell 
model has asB = 0.3910/0.4963, ess = 0.66944/5.726, 
and TUB = 15.035/77.106. The A - B interaction be- 
tween different molecules is determined by the Lorentz- 
Berthelot mixing rule. 

Rigid SPG/E water 59]: The water model is an isosceles 
triangle with sides roH — 1/3.166 and base hne 2roH * 
sin(109.47/2) (number of molecules studied: 1000). The 
oxygen-oxygen intermolecular interactions are given by 
the LJ pair potential (eoo = Ij foo = Ij and mo = 
15.9994/1.00794). There are no intermolecular LJ inter- 
actions for H-H or H-0. The three particles are charged 
with qo = -0.8476e/(47reo3.166Ao.650A:J/moZ)i/2 and 

1H = I90I/2 ensuring charge neutrality. 



Appendix B: How to delimit the first coordination 
shell? 



In all simulations the FCS cutoff was defined by plac- 
ing the cutoff at the first minimum of the RDF, which 
is the standard definition of the FCS for liquids [ij]. An 
alternative definition goes back to van der Waals [s^ . 
The FCS is here identified with a sphere of radius deter- 
mined by requiring that the average density within the 
FCS pintegrated equals the overall average density pmoan- 
Figure [201 shows for the KABLJ system the integrated 
local density of A particles calculated from the RDF (in- 
cluding the particle at the center) as a function of the 
distance to the origin. The "van der Waals distance" is 
slightly larger than the first minimum of the RDF. 
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FIG. 20. Integrated local density Pintegrated of the A parti- 
cles of the standard KABLJ mixture plotted as a function of 
distance to the origin (calculated from the A-particle RDF). 
The dashed black line gives the distance of the minimum of 
the RDF, the dashed magenta line gives the van der Waals 
distance, i.e., the distance at which the integrated local den- 
sity equals the overall average density of the system. The 
horizontal dashed blue line marks the mean density Pmcan of 
A particles in the simulated system {pmean = 0.8 * 1.2). 



FIG. 21. Effect on the A-particle ISF for the standard KABLJ 
mixture of different of ways of delimiting the FCS. (a) FCS 
is identified from the minimum beyond the first peak of the 
AA-particle RDF. (b) FCS is identified by the van der Waals 
distance, i.e., the distance at which the integrated local den- 
sity equals the mean density of the system. The van der Waals 
distance is slightly larger than the RDF minimum. Using the 
latter as defining the FCS cutoff radius gives a better repre- 
sentation of the dynamics, as seen from the inset. 



We applied this alternative definition of the FCS cutoff 
in Fig. [211 which shows the A-particle ISF for the (12, 6)- 
K ABL J mixture of Fig. simulated with, respectively, 
a cutoff at the first minimum of the RDF (Fig. [2lTa)) 
and the van der Waals cutoff (Fig. HHb)). 



Although the cutoff difference is merely 0.05, the van 
der Waals cutoff approximates better the reference ISF 
than does the RDF minimum cutoff. Thus it is possible 
that the van der Waals distance may serve as a better 
definition of the size of the FCS than the standard FCS 
definition. 




10" 10"' 10° lo' 10^ 10"^ 10* 



Identifying the size of the FCS for molecular systems is 
less straightforward, especially when different intermolec- 
ular interactions are involved. It is noteworthy how well 
the simple cutoff scheme in Fig. [TS] represents the dy- 
namics of the asymmetric dumbbell model. The slight 
deviations observed for the OTP model (Fig. [TWb)). 
disappear when the cutoff is increased from Vc — 1.47 
to Tc = 1.56 (Fig. [221) • This distance is close, but not 
identical, to the van der Waals distance calculated from 
the particle RDF (« 1.53). More work is needed to clar- 
ify the best way to delimit the FCS and define the FCS 
cutoff. 
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FIG. 22. Results for the ISF of the Lewis- Wahnstrom OTP 
model with a cutoff at rc = 1.56 (red), and a large reference 
cutoff (black). The deviations of Fig. [iTjb), in which — 
1.47, disappear by choosing this slightly larger cutoff, not far 
from the van der Waals cutoff (rc = 1.53). 
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